Entanglement Percolation in Quantum Complex Networks 
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Quantum networks are essential to quantum information distributed applications, and communi- 
cating over them is a key challenge. Complex networks have rich and intriguing properties, which 
are as yet unexplored in the quantum setting. Here, we study the effect of entanglement percolation 
as a means to establish long-distance entanglement between arbitrary nodes of quantum complex 
networks. We develop a theory to analytically study random graphs with arbitrary degree distri- 
bution and give exact results for some models. Our findings are in good agreement with numerical 
simulations and show that the proposed quantum strategies enhance the percolation threshold sub- 
stantially. Simulations also show a clear enhancement in small-world and other real-world networks. 
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A broad variety of natural and socioeconomic phenom- 
ena such as protein-protein interactions, the brain, the 
power-grid, friendship networks, the Internet and food- 
web, can all be modeled by graphs, i.e., by sets of nodes 
and edges representing the relation between them. Com- 
plex networks (CN) cover the wide range between regular 
lattices and completely random graphs. Their ubiquity 
has triggered an intense research activity on CN involving 
applied works, but also fundamental studies in theoreti- 
cal physics and mathematics that aim at unveiling their 
underlying principles. Understanding structural proper- 
ties of CN is very important as they crucially affect their 
functionality. For instance, the topology of a social net- 
work affects the spread of information or diseases, and 
the architecture of a computer network determines its 
robustness under router failures. The appearance of a 
giant connected component is a critical phenomenon re- 
lated to percolation theory P, Q ■ If, say, the infection 
probability of a disease exceeds a critical value, an out- 
burst of the disease occurs infecting most of the popu- 
lation, whereas below the critical value only a vanishing 
amount of nodes is affected. This threshold and the size 
of the outburst are key properties that strongly depend 
on the CN structure. 

This Letter merges the field of CN with the also 
new and rapidly developing field of quantum informa- 
tion (QI). Networks already have a prominent role in QI, 
e.g., in the context of measurement-based quantum com- 
putation Q , in the characterization of graph states [3] , 
and, more naturally, as the physical substrate of nodes 
and channels for multipartite quantum communication 
protocols. Such a "quantum internet" supports QI 
applications that fill the technological gap between the 
already available bipartite applications like quantum key 
distribution, and the appealing but still remote quan- 
tum computer. The reason for studying quantum CN, 
as opposed to regular lattices considered so far, is three- 
fold. First, it is very plausible that future quantum com- 
munication networks acquire a complex topology resem- 
bling that of existing networks. This can certainly be 



the case if methods are developed to use current commu- 
nication networks at a quantum regime, or if new quan- 
tum networks grow driven by similar sociopolitical mech- 
anisms. Second, CN have nontrivial topological features 
not found in regular lattices, e.g., robustness to random 
errors, that can be useful in the design of novel quantum 
networks. Third, from a more theoretical perspective, 
although few qubit states have been successfully charac- 
terized, an understanding of multipartite entanglement 
of many-particle systems is still lacking. As we will show 
here, some quantities which are exceedingly hard to cal- 
culate for states defined over regular lattices can be car- 
ried with ease in CN (which are defined through statisti- 
cal properties). 

The goal of entanglement percolation Q is to establish 
a maximally entangled state (which via quantum telepor- 
tation amounts to a perfect single-use quantum channel) 
between two arbitrary nodes of a network, where nodes 
are connected by partially entangled states. This can be 
easily achieved between two neighbor nodes since a par- 
tially entangled state can be converted into a maximally 
entangled state (for short singlet) with a probability that 
depends on the initial amount of entanglement. In ID 
networks, a singlet can only be established between two 
remote nodes with a probability that decays exponen- 
tially with the distance between them [6]. In contrast, 
higher dimensional networks have an increasing number 
of possible connecting paths giving rise to a percolation 
effect, i.e., to the appearance of a cluster spanning a sig- 
nificant fraction S of the network. Once a path of singlets 
is formed, a maximally entangled state can be established 
between the end-nodes by means of entanglement swap- 
ping. Thus, two arbitrary nodes can be connected by a 
singlet if both belong to the same cluster. For large net- 
works, this happens with a probability S'^ independent of 
the distance between the nodes and of the size of the net- 
work Ml, but that strongly depends on its topology. Aci'n 
et al. [61] found a type of lattice where the value of the 
percolation threshold, i.e., the required amount of initial 
entanglement, can be reduced by some quantum prepro- 
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cessing: local quantum measurements are performed on 
some selected nodes of the honeycomb lattice turning it 
into a triangular with a lower threshold. In other 
local strategies were studied in detail and entanglement 
percolation was shown to appear in few more 2D lattices. 

In this Letter, we study entanglement percolation in 
complex networks. We first derive general results for 
random networks with arbitrary uncorrelated degree dis- 
tributions. From here, we can obtain exact values for 
the percolation threshold and the probability S. We 
next complement these results with numerical simula- 
tions on these and other paradigmatic CN, with emphasis 
on those that mimic existing communication networks. 

A network is naturally represented by a graph, which 
consists of N vertices and a number of edges that connect 
some pairs of vertices. The degree of a vertex, fc, is the 
number of edges connecting to it, and in general follows a 
probability distribution P(fc). A connected component, 
or cluster, is a subgraph where any two vertices are con- 
nected by at least one path and no more vertices can be 
added without losing this property. If edges are occupied 
with probability p, one can ask which is the distribu- 
tion of connected component sizes when this probability 
varies. For low values of p many small clusters exist. As 
p increases clusters start to grow and join each other. In 
general the cluster size is vanishingly small compared to 
the size of the network, until a critical value oi p = Pc 
where a giant connected component (GCC) starts to span 
it: in the asymptotic limit, this is the only component of 
finite relative size 5* > for all p > Pc- 

We consider quantum networks, as in Ref. 0, whose 
edges correspond to pairs of pure entangled states 
where lip) = a/Aq |00) + \/A7|11) is a two qubit state 
with Schmidt coefficients \Ao > > 0. Each par- 

tially entangled state can be converted into a singlet 
by LOCC with the singlet conversion probability (SCP) 
p = min[l,2(l — Ao)] [Sj. If the two nodes share two 
copies of the probability that at least one of them 
is converted into a singlet is of course p2 — 2p — p^ . 
However, they can do better since the largest Schmidt 
cocfBcient of is Ao, and hence its optimal SCP is 

P2 = min[l, 2(1 — Aq)] = min(l, 2j)— p^/2). From here on, 
"with" and "without distillation" respectively mean that 
the optimal or the former sequential transformation is 
used. Applying the above singlet conversion to each edge 
of the network is called classical entanglement percola- 
tion 0. As already mentioned, a new quantum feature 
appears when we allow for other local quantum opera- 
tions that transform the network geometry and thereby 
change the percolation threshold. The basic ingredient is 
entanglement swapping (or "swap" ) illustrated in Fig. [1] 
The party at a central node c performs a Bell measure- 
ment on two qubits, each of them belonging to states 
shared with different nodes (a and b). After this opera- 
tion, the central qubits become disentangled from a and 
&, but in return an entangled state is established between 



a and 6, which on average has the same SCP as Note, 
however, that this operation cannot be repeated using 
newborn edges since they are not in state The "g- 
swap" operation that we use here performs the "swap" 
transformation between successive pairs of neighbors of 
a central node of degree g, effectively changing the initial 
g-star into a g-cycle (Fig. [1]). These g-swaps are made 
on nodes of certain degrees, depending on the network 
structure with the only condition that no operation is 
done on neighboring sites. Thus, the strategy uses only 
local information of the network: the degree of the target 
node and the status of its neighbors. 

We now proceed to present our results for general ran- 
dom graphs of arbitrary (uncorrelated) degree distribu- 
tion. We use the generating function formalism l^, 11 1. 
We begin by defining the generating function of the de- 
gree distribution by Go[x) — Y^^=aPi.^)^^ ■ Since P{k) 
is a probability. Go is well normalized, Go(l) = 1, and 
convergent for \x\ < 1. What makes generating functions 
specially suited here is that it is straightforward to convo- 
lute probability distributions; e.g., Go(a;)'" generates the 
probability that m arbitrary nodes have a total of k edges 
leading to them. The key probability distribution in per- 
colation is that of reaching a finite component of size s 
when following a random edge to one of its ends, gener- 
ated by Hi (x) . The size is zero when the edge is not occu- 
pied. With probability p2 it is occupied, and the discov- 
ered vertex has degree distribution Pi(fc) — kP{k)/{k) 

generated by X]fc°=o ^i(^)-^'^ ~ ^ = xGi{x). Thus, 
s is one plus the size of the new fc — 1 reachable com- 
ponents. The size of each of them is again generated by 
iJi(x), and their sum by E^o There- 
fore, Hi{x) fulfills the consistency equation 



Hiix) = {l^p2)+P2xGiiHi{x)). 



(1) 



It is crucial to notice that by restricting to finite s 
we have explicitly excluded the infinite GCC from Hi (x) 
[l^ . The total probability that an edge connects to a 
finite component is u = Hi{l), and due to ^ satisfies 
u — 1 — p2 + p2Gi{u). Below the critical p < pc there is 
a unique solution u = I, while at p > p^ & new solution 
u < 1 appears. Hence, the critical probability is the 
smallest value of p for which the second solution appears. 

We can similarly compute the function Hq (x) generat- 
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FIG. 1: Entanglement swapping (left) and 5-swap transforma- 
tion (right). Dots represent qubits, circles nodes, solid lines 
partially entangled states \tp), and dashed lines correspond to 
states with the same SCP as \tp}. 
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ing the probability that a randomly chosen vertex be- 
longs to a finite component of size s. Since k edges 
emerge from this vertex with probability P{k) we find 
Ho{x) — xGq[Hi(x)]. The probability that a randomly 
chosen vertex belongs to the GCC is 5 = 1 — Ho{l) ~ 
1 — Go{u). From Ho{x), we can compute many other rel- 
evant properties, e.g., the mean component size below pc, 
(s) = Hq{1), which behaves like a susceptibility diverging 
at the phase transition. 

Our goal now is to understand how quantum opera- 
tions change the percolation properties of the original 
network. Every particular g-swap can be implemented 
(or not) with probability 11^ (or l—Hq) on nodes of degree 
q. Giving the values for each Ilq specifies the quantum 
strategy. Now, instead of arriving to a vertex of degree q 
connecting to other q — 1 components, after a g-swap op- 
eration we arrive to a cycle of q nodes (including the one 
we are coming from) connected via links with SCP equal 
to p [l^. The size of such cycle and its emerging com- 
ponents is the sum of the probabilities that a connected 
string of I vertices in the cycle is reached, multiplied by 
the generating function of the emerging component sizes: 

C,{x) - ^(/ + l)p'(l~p)2 xGi(i?i(x)) 



1=0 
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+ [qp''-Hl [xGi{H,{x))\ .(2) 

Therefore, the new Hi (x) is of the same form of Eq. ^ 
plus a term Hi^q{x) for each g-swap: 

Hi{x) = l-p2+P2xGiiHi{x)) + J2^qHi^gix) (3) 

q>2 



i?i,,(x) = Pi{q) 



{P2 - 1) -P2X {Hi{x) 
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Cq{x) 



At this stage we can already calculate pc as the value of 
p for which there exists a solution < u = Hi{l) < 1 
to ([3]) at X = 1. It is easy to convince oneself that each 
separate contribution Hi^q{l) either increases or lowers 
the percolation threshold and therefore for the optimal 
strategy each Hq is either or 1. 

In order to calculate the relative jize of the giant com- 
ponent S we need to find the new Hq(x). Given a vertex 
q, there is a probability rjq that a g-swap can be per- 
formed on it, changing its degree from q to zero and hence 



Hq(x) = xG[ 



{Hi{x))^ 



q>2 



1 



(i/i(x))' 



From here, the order parameter 5 = 1 — Ho{l) can be 
readily obtained by plugging in the solution of Hi(l) = u. 
The probability rjq that no other q'-swap was performed 
in its neighbors can be computed up to any order in 
p exactly [l2|. When only one type of g-swap is per- 
formed, a compact, approximate expression exists rjq « 



(1 -I- |Pi(g))~^. In the following, we study entanglement 
percolation enhancement in particular complex network 
models. 

The Bethe lattice is an infinite regular tree where every 
vertex has the same degree. Inserting Go{x) = x'^ and 
Gi{x) — x'^^^ into 111 and Q and solving the point at 
which the solution u < 1 appears gives (q — 1)^^ ^ P2 
and (<? - l)-i = (1 - p)-^{2p + p<i[p{q - 1) - (g + 1)]} 
after g-swap is applied. Therefore, g-swap gives always a 
better threshold except for an infinite ID chain (g = 2). 

Erdos-Renyi (ER) graphs are maximally random 
graphs with the only constraint (fc) ~ z. An ER net- 
work has N vertices and each pair of them holds an 
edge with probability z/N. The degrees follow a Pois- 
son distribution and Go{x) = Gi{x) = exp[z(a; — 1)]. 
The original threshold is given by p2 = 1/z. After, 
e.g., the 2-swap and 3-swap operations, the thresholds 
are respectively j=P2 + P2 + z{2p — p'^)] and 

^ = P2 + ze~^[—p2 + z{l + p — p"^)]. In general the 
performance of different g-swaps depends on the mean 
degree z. Figure [5^ shows the gain 7 — {pc —pc)/Pc- 
The threshold can improve by a 20% (24% without dis- 
tillation), compared to the 3% (16% without distillation) 
of the honeycomb-triangular transformation in The 
size of the GCC can be calculated exactly from the above 
results [l2| and is plotted in Fig. [^^b) showing perfect 
agreement with numerical simulations. Note that S at- 
tains smaller values for the transformed network even at 
p = 1 due to the target nodes that become detached after 
the g-swap. Since one can always choose where to start 
applying the g-swaps, the probability of connecting two 
arbitrary nodes is actually , where S = SSi / Si is the 
probability that a node belongs to the GCC excluding 
the detached nodes and Si is the GCC size at p = 1. 
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FIG. 2; (a) Gain 7 in the ER model as a function of the mean 
degree z after 2-swap (dotted), 2, 3-swap (dashed), 2,3,4-swap 
(dot-dashed) and optimal g-swaps (solid), (b) Size of the GCC 
as a function of p in the z = 2.5 ER network without distilla- 
tion before (squares) and after 2, 3-swap (circles). Points are 
simulation results for N = 10^. (c) Same as (a) for scale- free 
model with r = 1 as a function of k. (d) Same as (b) for 
WWW [il with cutoff fc = 15 (AT ~ 10^). 
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FIG. 3: Numerical results for the SW model with 2-swap and 
distillation, (a) pc as a function of (p before (squares) and 
after (disks). Solid line is the classical analytic solution from 
[m . (b) Size of the GCC as a function of p before (squares) 
and after (circles) operations for (p — 0.25. 



Scale-free networks: Real world networks are not Pois- 
sonian but typically have a power-law (scale-free) degree 
distribution. Such networks have a GCC at all p, and a 
threshold appears only in scenarios where there is a cut- 
off at high degrees due to, e.g., targeted attacks, physical 
constraints or saturation effects. We have taken an expo- 
nential cutoff, P{k) — C/c^'^e"'^/", and find once again 
a significant improvement in pc, up to 25% with distil- 
lation [Fig. [2tc)]. We have performed simulations for 
a real- world WWW network [13|, where we have intro- 
duced a cutoff neglecting nodes with A: > 15 and, de- 
spite the finite-size effects, we clearly see a decrease in pc 
[Fig.Hd)]. 

A Small World (SW) flil is a network with ordered 
local structure and high level of clustering but still with 
surprisingly low average path length. We consider an 
SW model generated by placing N vertices in a ID chain. 
Then N additional random edges, "shortcuts" , are added 
with probability (f>. In this case, degrees are strongly cor- 
related and thus the above analysis does not apply. Still, 
simulations show that g-swaps can decrease Pc (Fig. [3]) . 

To summarize, we have studied entanglement perco- 
lation in complex networks — which are defined through 
their statistical properties and whose exact structure is 
not necessarily known — and shown that by a strategy 
that only uses local information we can decrease the per- 
colation threshold substantially (in some cases almost 1 
order of magnitude more than found on the honeycomb 
lattice). We have exactly solved models of random graphs 
with arbitrary uncorrelated degree distribution. More- 
over, we have numerically explored the behavior in SW 
networks and the WWW. In all cases, we find that entan- 
glement percolation, far from being restricted to partic- 
ular networks, is a quite general feature. Entanglement 
percolation is interesting from the discipline of CN since 
quantum mechanics introduces a new paradigm where 
networks can suffer non-trivial structural changes before 
the percolation process starts. We believe our results 
open prospects for many interesting synergies between 
the fields of quantum information and complex networks. 

An essential issue for entanglement percolation 
schemes is their stability under noise. For amplitude 



damping, one can adapt the prescription in 15|. More 
general scenarios might require schemes that, e.g., use 
multipartite entanglement, use global knowledge of the 
network, concentrate different paths in a single one, or 
allow for the generation of new (imperfect) bonds during 
the protocol. Finally, note that CN can enjoy much bet- 
ter stability than lattices. In particular, in many CN, the 
average path length scales with the size of the network 
as log N, as opposed to N^^"^ for d-dimensional lattices. 
Thus, the probability of establishing entanglement be- 
tween two arbitrary nodes has a linear falloff, as opposed 
to a fatal exponential decay, and hence is non-negligible 
in large, but finite, networks. 

We are grateful to M. Boguna and A.V. Goltsev for 
discussions. We acknowledge financial support from the 
Spanish MICINN, through the Ramon y Cajal (J. C.) and 
AP2008-03048 (M. C.) and projects FIS2008-01236 and 
QOIT (CONSOLIDER2006-00019), and the Generahtat 
de Catalunya Contract No. 2009SGR985. 

Note added. — Soon after the completion of this work, 
in 16] Perseguers et al. showed radically different proper- 
ties between quantum and classical Erdos-Renyi graphs, 
and also stressed the interest in quantum CN. Interesting 
results for entanglement distribution in regular lattices 
have also been recently reported in (l7j . 
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